Genomic Profiling to Contextualize the Results of Intervention for Smoldering Multiple Myeloma.

PURPOSE: Early intervention for High-Risk Smoldering Multiple Myeloma (HR-SMM) achieves deep and prolonged responses. It is unclear if beneficial outcomes are due to treatment of less complex, susceptible disease or inaccuracy in clinical definition of cases entered. EXPERIMENTAL DESIGN: Here, we interrogated whole genome and whole exome sequencing for 54 patients across two HR-SMM interventional studies (NCT01572480, NCT02279394). RESULTS: We reveal that the genomic landscape of treated HR-SMM is generally simple as compared to Newly Diagnosed (ND)MM counterparts with less inactivation of tumor suppressor genes, RAS pathway mutations, MYC disruption, and APOBEC contribution. The absence of these events parallels that of indolent precursor conditions, possibly explaining overall excellent outcomes. However, some patients harboring genomic complexity fail to sustain response and experience resistant, progressive disease. Overall, clinical risk scores do not effectively discriminate between genomically indolent and aggressive disease. CONCLUSIONS: Genomic profiling can contextualize the advantage of early intervention in SMM and guide personalization of therapy.

Whole-exome sequencing is feasible on a fresh-frozen skin sample of intravascular large B cell lymphoma.

Intravascular large B-cell lymphoma (IVLBCL) is a rare aggressive extranodal non-Hodgkin lymphoma. The predominant, if not exclusive, growth of neoplastic cells within the lumina of small-sized vessels represents the hallmark of the disease. Diagnosis is challenging due to the absence of marked lymphadenopathy, the highly heterogeneous clinical presentation, and the rarity of the condition. Clinical presentation is characterized by variable combinations of nonspecific signs and symptoms (such as fever and weight loss), organ-specific focal manifestations due to altered perfusion, and hemophagocytic syndrome. The rarity of this entity and the paucity of neoplastic cells in biopsy samples hamper the study of recurrent molecular abnormalities. The purpose of this study was to explore the feasibility of a different approach to recover a sufficient amount of DNA of acceptable quality to perform next-generation sequencing studies. Here, we report the findings of whole-exome next-generation sequencing performed on a fresh-frozen cutaneous sample of IVLBCL, paired with the patient saliva used as germline DNA. To increase the cancer cell fraction, only the subcutaneous tissue was selected. With this approach, we obtained high-quality DNA and were able to identify oncogenic mutations specific for this entity and recapitulating its post-germinal center origin, even if the tumor fraction was low. Molecular studies performed on fresh-frozen cutaneous sample are feasible in IVLBCL, especially when analysis is restricted to the subcutaneous tissue. Wide adoption of this reproducible and cost-effective approach may foster further studies, which may be of help in supporting diagnosis, providing pathogenetic insights, and guiding treatment decisions.

Genomic Classification and Individualized Prognosis in Multiple Myeloma.

PURPOSE: Outcomes for patients with newly diagnosed multiple myeloma (NDMM) are heterogenous, with overall survival (OS) ranging from months to over 10 years. METHODS: To decipher and predict the molecular and clinical heterogeneity of NDMM, we assembled a series of 1,933 patients with available clinical, genomic, and therapeutic data. RESULTS: Leveraging a comprehensive catalog of genomic drivers, we identified 12 groups, expanding on previous gene expression-based molecular classifications. To build a model predicting individualized risk in NDMM (IRMMa), we integrated clinical, genomic, and treatment variables. To correct for time-dependent variables, including high-dose melphalan followed by autologous stem-cell transplantation (HDM-ASCT), and maintenance therapy, a multi-state model was designed. The IRMMa model accuracy was significantly higher than all comparator prognostic models, with a c-index for OS of 0.726, compared with International Staging System (ISS; 0.61), revised-ISS (0.572), and R2-ISS (0.625). Integral to model accuracy was 20 genomic features, including 1q21 gain/amp, del 1p, TP53 loss, NSD2 translocations, APOBEC mutational signatures, and copy-number signatures (reflecting the complex structural variant chromothripsis). IRMMa accuracy and superiority compared with other prognostic models were validated on 256 patients enrolled in the GMMG-HD6 (ClinicalTrials.gov identifier: NCT02495922) clinical trial. Individualized patient risks were significantly affected across the 12 genomic groups by different treatment strategies (ie, treatment variance), which was used to identify patients for whom HDM-ASCT is particularly effective versus patients for whom the impact is limited. CONCLUSION: Integrating clinical, demographic, genomic, and therapeutic data, to our knowledge, we have developed the first individualized risk-prediction model enabling personally tailored therapeutic decisions for patients with NDMM.

Impact of Rare Structural Variant Events in Newly Diagnosed Multiple Myeloma.

PURPOSE: Whole-genome sequencing (WGS) of patients with newly diagnosed multiple myeloma (NDMM) has shown recurrent structural variant (SV) involvement in distinct regions of the genome (i.e., hotspots) and causing recurrent copy-number alterations. Together with canonical immunoglobulin translocations, these SVs are recognized as "recurrent SVs." More than half of SVs were not involved in recurrent events. The significance of these "rare SVs" has not been previously examined. EXPERIMENTAL DESIGN: In this study, we utilize 752 WGS and 591 RNA sequencing data from patients with NDMM to determine the role of rare SVs in myeloma pathogenesis. RESULTS: Ninety-four percent of patients harbored at least one rare SV event. Rare SVs showed an SV class-specific enrichment within genes and superenhancers associated with outlier gene expression. Furthermore, known myeloma driver genes recurrently impacted by point mutations were dysregulated by rare SVs. CONCLUSIONS: Overall, we demonstrate the association of rare SVs with aberrant gene expression supporting a potential driver role in myeloma pathogenesis.

Genomic and immune determinants of resistance to anti-CD38 monoclonal antibody-based therapy in relapsed refractory multiple myeloma.

Anti-CD38 antibody therapies have transformed multiple myeloma (MM) treatment. However, a large fraction of patients inevitably relapses. To understand this, we investigated 32 relapsed MM patients treated with daratumumab, lenalidomide, and dexamethasone (Dara-Rd; NCT03848676 ). Whole genome sequencing (WGS) before and after treatment pinpointed genomic drivers associated with early progression, including RPL5 loss and APOBEC mutagenesis. Flow cytometry on 202 blood samples, collected every three months until progression for 31 patients, revealed distinct immune changes significantly impacting clinical outcomes. Progressing patients exhibited significant depletion of CD38+ NK cells, persistence of T cell exhaustion, and reduced depletion of T-reg cells over time. These findings underscore the influence of immune composition and daratumumab-induced immune changes in promoting MM resistance. Integrating genomics and flow cytometry unveiled associations between adverse genomic features and immune patterns. Overall, this study sheds light on the intricate interplay between genomic complexity and the immune microenvironment driving resistance to Dara-Rd.

Genomic and immune signatures predict clinical outcome in newly diagnosed multiple myeloma treated with immunotherapy regimens.

Despite improving outcomes, 40% of patients with newly diagnosed multiple myeloma treated with regimens containing daratumumab, a CD38-targeted monoclonal antibody, progress prematurely. By integrating tumor whole-genome and microenvironment single-cell RNA sequencing from upfront phase 2 trials using carfilzomib, lenalidomide and dexamethasone with daratumumab ( NCT03290950 ), we show how distinct genomic drivers including high APOBEC mutational activity, IKZF3 and RPL5 deletions and 8q gain affect clinical outcomes. Furthermore, evaluation of paired bone marrow profiles, taken before and after eight cycles of carfilzomib, lenalidomide and dexamethasone with daratumumab, shows that numbers of natural killer cells before treatment, high T cell receptor diversity before treatment, the disappearance of sustained immune activation (that is, B cells and T cells) and monocyte expansion over time are all predictive of sustained minimal residual disease negativity. Overall, this study provides strong evidence of a complex interplay between tumor cells and the immune microenvironment that is predictive of clinical outcome and depth of treatment response in patients with newly diagnosed multiple myeloma treated with highly effective combinations containing anti-CD38 antibodies.

The Vk*MYC Mouse Model recapitulates human multiple myeloma evolution and genomic diversity.

Despite advancements in profiling multiple myeloma (MM) and its precursor conditions, there is limited information on mechanisms underlying disease progression. Clincal efforts designed to deconvolute such mechanisms are challenged by the long lead time between monoclonal gammopathy and its transformation to MM. MM mouse models represent an opportunity to overcome this temporal limitation. Here, we profile the genomic landscape of 118 genetically engineered Vk*MYC MM and reveal that it recapitulates the genomic heterogenenity and life history of human MM. We observed recurrent copy number alterations, structural variations, chromothripsis, driver mutations, APOBEC mutational activity, and a progressive decrease in immunoglobulin transcription that inversely correlates with proliferation. Moreover, we identified frequent insertional mutagenesis by endogenous retro-elements as a murine specific mechanism to activate NF-kB and IL6 signaling pathways shared with human MM. Despite the increased genomic complexity associated with progression, advanced tumors remain dependent on MYC expression, that drives the progression of monoclonal gammopathy to MM.

Mechanisms of antigen escape from BCMA- or GPRC5D-targeted immunotherapies in multiple myeloma.

B cell maturation antigen (BCMA) target loss is considered to be a rare event that mediates multiple myeloma (MM) resistance to anti-BCMA chimeric antigen receptor T cell (CAR T) or bispecific T cell engager (TCE) therapies. Emerging data report that downregulation of G-protein-coupled receptor family C group 5 member D (GPRC5D) protein often occurs at relapse after anti-GPRC5D CAR T therapy. To examine the tumor-intrinsic factors that promote MM antigen escape, we performed combined bulk and single-cell whole-genome sequencing and copy number variation analysis of 30 patients treated with anti-BCMA and/or anti-GPRC5D CAR T/TCE therapy. In two cases, MM relapse post-TCE/CAR T therapy was driven by BCMA-negative clones harboring focal biallelic deletions at the TNFRSF17 locus at relapse or by selective expansion of pre-existing subclones with biallelic TNFRSF17 loss. In another five cases of relapse, newly detected, nontruncating, missense mutations or in-frame deletions in the extracellular domain of BCMA negated the efficacies of anti-BCMA TCE therapies, despite detectable surface BCMA protein expression. In the present study, we also report four cases of MM relapse with biallelic mutations of GPRC5D after anti-GPRC5D TCE therapy, including two cases with convergent evolution where multiple subclones lost GPRC5D through somatic events. Immunoselection of BCMA- or GPRC5D-negative or mutant clones is an important tumor-intrinsic driver of relapse post-targeted therapies. Mutational events on BCMA confer distinct sensitivities toward different anti-BCMA therapies, underscoring the importance of considering the tumor antigen landscape for optimal design and selection of targeted immunotherapies in MM.

Molecular Evolution of Classic Hodgkin Lymphoma Revealed Through Whole-Genome Sequencing of Hodgkin and Reed Sternberg Cells.

The rarity of malignant Hodgkin and Reed Sternberg (HRS) cells in classic Hodgkin lymphoma (cHL) limits the ability to study the genomics of cHL. To circumvent this, our group has previously optimized fluorescence-activated cell sorting to purify HRS cells. Using this approach, we now report the whole-genome sequencing landscape of HRS cells and reconstruct the chronology and likely etiology of pathogenic events leading to cHL. We identified alterations in driver genes not previously described in cHL, APOBEC mutational activity, and the presence of complex structural variants including chromothripsis. We found that high ploidy in cHL is often acquired through multiple, independent chromosomal gains events including whole-genome duplication. Evolutionary timing analyses revealed that structural variants enriched for RAG motifs, driver mutations in B2M, BCL7A, GNA13, and PTPN1, and the onset of AID-driven mutagenesis usually preceded large chromosomal gains. This study provides a temporal reconstruction of cHL pathogenesis. SIGNIFICANCE: Previous studies in cHL were limited to coding sequences and therefore not able to comprehensively decipher the tumor complexity. Here, leveraging cHL whole-genome characterization, we identify driver events and reconstruct the tumor evolution, finding that structural variants, driver mutations, and AID mutagenesis precede chromosomal gains. This article is highlighted in the In This Issue feature, p. 171.

Tracking the evolution of therapy-related myeloid neoplasms using chemotherapy signatures.

Patients treated with cytotoxic therapies, including autologous stem cell transplantation, are at risk for developing therapy-related myeloid neoplasms (tMN). Preleukemic clones (ie, clonal hematopoiesis [CH]) are detectable years before the development of these aggressive malignancies, although the genomic events leading to transformation and expansion are not well defined. Here, by leveraging distinctive chemotherapy-associated mutational signatures from whole-genome sequencing data and targeted sequencing of prechemotherapy samples, we reconstructed the evolutionary life-history of 39 therapy-related myeloid malignancies. A dichotomy was revealed, in which neoplasms with evidence of chemotherapy-induced mutagenesis from platinum and melphalan were hypermutated and enriched for complex structural variants (ie, chromothripsis), whereas neoplasms with nonmutagenic chemotherapy exposures were genomically similar to de novo acute myeloid leukemia. Using chemotherapy-associated mutational signatures as temporal barcodes linked to discrete clinical exposure in each patient's life, we estimated that several complex events and genomic drivers were acquired after chemotherapy was administered. For patients with prior multiple myeloma who were treated with high-dose melphalan and autologous stem cell transplantation, we demonstrate that tMN can develop from either a reinfused CH clone that escapes melphalan exposure and is selected after reinfusion, or from TP53-mutant CH that survives direct myeloablative conditioning and acquires melphalan-induced DNA damage. Overall, we revealed a novel mode of tMN progression that is not reliant on direct mutagenesis or even exposure to chemotherapy. Conversely, for tMN that evolve under the influence of chemotherapy-induced mutagenesis, distinct chemotherapies not only select preexisting CH but also promote the acquisition of recurrent genomic drivers.

Impact of rare structural variant events in newly diagnosed multiple myeloma.

Whole genome sequencing (WGS) of newly diagnosed multiple myeloma patients (NDMM) has shown recurrent structural variant (SV) involvement in distinct regions of the genome (i.e. hotspots) and causing recurrent copy number alterations. Together with canonical immunoglobulin translocations, these SVs are recognized as "recurrent SVs". More than half SVs were not involved in recurrent events. The significance of these "rare SVs" has not been previously examined. In this study, we utilize 752 WGS and 591 RNA-seq data from NDMM patients to determine the role of rare SVs in myeloma pathogenesis. 94% of patients harbored at least one rare SV event. Rare SVs showed an SV-class specific enrichment within genes and superenhancers associated with outlier gene expression. Furthermore, known myeloma driver genes recurrently impacted by point mutations were dysregulated by rare SVs. Overall, we demonstrate the association of rare SVs with aberrant gene expression supporting a driver role in myeloma pathogenesis. SIGNIFICANCE: Characterization of multiple myeloma genome revealed that more than half structural variants are not involved in recurrent events. Here, we demonstrate that these rare SVs hold potential for myeloma pathogenesis through their gene expression impact. Rare SVs contribute to MM heterogeneity and have implications for development of individualized treatment.

Waldenström macroglobulinemia whole genome reveals prolonged germinal center activity and late copy number aberrations.

The genomic landscape of Waldenström macroglobulinemia (WM) is characterized by somatic mutations in MYD88, present from the precursor stages. Using the comprehensive resolution of whole genome sequencing (WGS) in 14 CD19-selected primary WM samples; comparing clonal and subclonal mutations revealed that germinal center (GC) mutational signatures SBS9 (poly-eta) and SBS84 (AID) have sustained activity, suggesting that the interaction between WM and the GC continues over time. Expanding our cohort size with 33 targeted sequencing samples, we interrogated the WM copy number aberration (CNA) landscape and chronology. Of interest, CNA prevalence progressively increased in symptomatic WM and relapsed disease when compared with stable precursor stages, with stable precursors lacking genomic complexity. Two MYD88 wild-type WGS contained a clonal gain affecting chromosome 12, which is typically an early event in chronic lymphocytic leukemia. Molecular time analysis demonstrated that both chromosomal 12 gain events occurred early in cancer development whereas other CNA changes tend to occur later in the disease course and are often subclonal. In summary, WGS analysis in WM allows the demonstration of sustained GC activity over time and allows the reconstruction of the temporal evolution of specific genomic features. In addition, our data suggest that, although MYD88-mutations are central to WM clone establishment and can be observed in precursor disease, CNA may contribute to later phases, and may be used as a biomarker for progression risk from precursor conditions to symptomatic disease.

Whole-genome sequencing reveals complex genomic features underlying anti-CD19 CAR T-cell treatment failures in lymphoma.

CD19-directed chimeric antigen receptor (CAR-19) T cells are groundbreaking immunotherapies approved for use against large B-cell lymphomas. Although host inflammatory and tumor microenvironmental markers associate with efficacy and resistance, the tumor-intrinsic alterations underlying these phenomena remain undefined. CD19 mutations associate with resistance but are uncommon, and most patients with relapsed disease retain expression of the wild-type receptor, implicating other genomic mechanisms. We therefore leveraged the comprehensive resolution of whole-genome sequencing to assess 51 tumor samples from 49 patients with CAR-19-treated large B-cell lymphoma. We found that the pretreatment presence of complex structural variants, APOBEC mutational signatures, and genomic damage from reactive oxygen species predict CAR-19 resistance. In addition, the recurrent 3p21.31 chromosomal deletion containing the RHOA tumor suppressor was strongly enriched in patients for whom CAR T-cell therapy failed. Pretreatment reduced expression or monoallelic loss of CD19 did not affect responses, suggesting CAR-19 therapy success and resistance are related to multiple mechanisms. Our study showed that tumor-intrinsic genomic alterations are key among the complex interplay of factors that underlie CAR-19 efficacy and resistance for large B-cell lymphomas.

Functional Impact of Genomic Complexity on the Transcriptome of Multiple Myeloma.

PURPOSE: Multiple myeloma is a biologically heterogenous plasma-cell disorder. In this study, we aimed at dissecting the functional impact on transcriptome of gene mutations, copy-number abnormalities (CNA), and chromosomal rearrangements (CR). Moreover, we applied a geno-transcriptomic approach to identify specific biomarkers for personalized treatments. EXPERIMENTAL DESIGN: We analyzed 514 newly diagnosed patients from the IA12 release of the CoMMpass study, accounting for mutations in multiple myeloma driver genes, structural variants, copy-number segments, and raw-transcript counts. We performed an in silico drug sensitivity screen (DSS), interrogating the Cancer Dependency Map (DepMap) dataset after anchoring cell lines to primary tumor samples using the Celligner algorithm. RESULTS: Immunoglobulin translocations, hyperdiploidy and chr(1q)gain/amps were associated with the highest number of deregulated genes. Other CNAs and specific gene mutations had a lower but very distinct impact affecting specific pathways. Many recurrent genes showed a hotspot (HS)-specific effect. The clinical relevance of double-hit multiple myeloma found strong biological bases in our analysis. Biallelic deletions of tumor suppressors and chr(1q)-amplifications showed the greatest impact on gene expression, deregulating pathways related to cell cycle, proliferation, and expression of immunotherapy targets. Moreover, our in silico DSS showed that not only t(11;14) but also chr(1q)gain/amps and CYLD inactivation predicted differential expression of transcripts of the BCL2 axis and response to venetoclax. CONCLUSIONS: The multiple myeloma genomic architecture and transcriptome have a strict connection, led by CNAs and CRs. Gene mutations impacted especially with HS-mutations of oncogenes and biallelic tumor suppressor gene inactivation. Finally, a comprehensive geno-transcriptomic analysis allows the identification of specific deregulated pathways and candidate biomarkers for personalized treatments in multiple myeloma.

Copy number signatures predict chromothripsis and clinical outcomes in newly diagnosed multiple myeloma.

Chromothripsis is detectable in 20-30% of newly diagnosed multiple myeloma (NDMM) patients and is emerging as a new independent adverse prognostic factor. In this study we interrogate 752 NDMM patients using whole genome sequencing (WGS) to investigate the relationship of copy number (CN) signatures to chromothripsis and show they are highly associated. CN signatures are highly predictive of the presence of chromothripsis (AUC = 0.90) and can be used identify its adverse prognostic impact. The ability of CN signatures to predict the presence of chromothripsis is confirmed in a validation series of WGS comprised of 235 hematological cancers (AUC = 0.97) and an independent series of 34 NDMM (AUC = 0.87). We show that CN signatures can also be derived from whole exome data (WES) and using 677 cases from the same series of NDMM, we are able to predict both the presence of chromothripsis (AUC = 0.82) and its adverse prognostic impact. CN signatures constitute a flexible tool to identify the presence of chromothripsis and is applicable to WES and WGS data.

Specific targeting of the KRAS mutational landscape in myeloma as a tool to unveil the elicited antitumor activity.

Alterations in KRAS have been identified as the most recurring somatic variants in the multiple myeloma (MM) mutational landscape. Combining DNA and RNA sequencing, we studied 756 patients and observed KRAS as the most frequently mutated gene in patients at diagnosis; in addition, we demonstrated the persistence or de novo occurrence of the KRAS aberration at disease relapse. Small-molecule inhibitors targeting KRAS have been developed; however, they are selective for tumors carrying the KRASG12C mutation. Therefore, there is still a need to develop novel therapeutic approaches to target the KRAS mutational events found in other tumor types, including MM. We used AZD4785, a potent and selective antisense oligonucleotide that selectively targets and downregulates all KRAS isoforms, as a tool to dissect the functional sequelae secondary to KRAS silencing in MM within the context of the bone marrow niche and demonstrated its ability to significantly silence KRAS, leading to inhibition of MM tumor growth, both in vitro and in vivo, and confirming KRAS as a driver and therapeutic target in MM.

mmsig: a fitting approach to accurately identify somatic mutational signatures in hematological malignancies.

Mutational signatures have emerged as powerful biomarkers in cancer patients, with prognostic and therapeutic implications. Wider clinical utility requires access to reproducible algorithms, which allow characterization of mutational signatures in a given tumor sample. Here, we show how mutational signature fitting can be applied to hematological cancer genomes to identify biologically and clinically important mutational processes, highlighting the importance of careful interpretation in light of biological knowledge. Our newly released R package mmsig comes with a dynamic error-suppression procedure that improves specificity in important clinical and biological settings. In particular, mmsig allows accurate detection of mutational signatures with low abundance, such as those introduced by APOBEC cytidine deaminases. This is particularly important in the most recent mutational signature reference (COSMIC v3.1) where each signature is more clearly defined. Our mutational signature fitting algorithm mmsig is a robust tool that can be implemented immediately in the clinic.

A Journey Through Myeloma Evolution: From the Normal Plasma Cell to Disease Complexity.

The knowledge of cancer origin and the subsequent tracking of disease evolution represent unmet needs that will soon be within clinical reach. This will provide the opportunity to improve patient's stratification and to personalize treatments based on cancer biology along its life history. In this review, we focus on the molecular pathogenesis of multiple myeloma (MM), a hematologic malignancy with a well-known multi-stage disease course, where such approach can sooner translate into a clinical benefit. We describe novel insights into modes and timing of disease initiation. We dissect the biology of the preclinical and pre-malignant phases, elucidating how knowledge of the genomics of the disease and the composition of the microenvironment allow stratification of patients based on risk of disease progression. Then, we explore cell-intrinsic and cell-extrinsic drivers of MM evolution to symptomatic disease. Finally, we discuss how this may relate to the development of refractory disease after treatment. By integrating an evolutionary view of myeloma biology with the recent acquisitions on its clonal heterogeneity, we envision a way to drive the clinical management of the disease based on its detailed biological features more than surrogates of disease burden.

Early Relapse Risk in Patients with Newly Diagnosed Multiple Myeloma Characterized by Next-generation Sequencing.

PURPOSE: Duration of first remission is important for the survival of patients with multiple myeloma. EXPERIMENTAL DESIGN: From the CoMMpass study (NCT01454297), 926 patients with newly diagnosed multiple myeloma, characterized by next-generation sequencing, were analyzed to evaluate those who experienced early progressive disease (PD; time to progression, TTP ≤18 months). RESULTS: After a median follow-up of 39 months, early PD was detected in 191/926 (20.6%) patients, 228/926 (24.6%) patients had late PD (TTP >18 months), while 507/926 (54.8%) did not have PD at the current follow-up. Compared with patients with late PD, patients with early PD had a lower at least very good partial response rate (47% vs. 82%, P < 0.001) and more frequently acquired double refractoriness to immunomodulatory drugs (IMiD) and proteasome inhibitors (PI; 21% vs. 8%, P < 0.001). Patients with early PD were at higher risk of death compared with patients with late PD and no PD (HR, 3.65; 95% CI, 2.7-4.93; P < 0.001), showing a dismal median overall survival (32.8 months). In a multivariate logistic regression model, independent factors increasing the early PD risk were TP53 mutation (OR, 3.78, P < 0.001), high lactate dehydrogenase levels (OR, 3.15, P = 0.006), λ-chain translocation (OR, 2.25, P = 0.033), and IGLL5 mutation (OR, 2.15, P = 0.007). Carfilzomib-based induction (OR, 0.15, P = 0.014), autologous stem-cell transplantation (OR, 0.27, P < 0.001), and continuous therapy with PIs and IMiDs (OR, 0.34, P = 0.024) mitigated the risk of early PD. CONCLUSIONS: Early PD identifies a high-risk multiple myeloma population. Further research is needed to better identify baseline features predicting early PD and the optimal treatment approaches for patients at risk.